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We investigate spin and charge current through a quantum dot pumped by a time- varying mag- 
netic field. Using the density matrix method, quantum rate equations for the electronic occupation 
numbers in the quantum dot are obtained and solved in the stationary state limit for a wide set of 
setup parameters. Both charge and spin current are expressed explicitly in terms of several relevant 
parameters and analyzed in detail. The results suggest a way of optimizing experimental setup 
parameters to obtain an maximal spin current without the charge current flow. 

PACS numbers: 



Spintronics is an emerging active research field, which is based on the effective control of electron spin in addition 
to its charge degree of freedom^. One common operational principle for spin-based devices is how to generate a spin 
current. It is known that a spin current is usually accompanied by a charge current, which generates heat. And this 
may pose a severe problem in microelectronics as the devices become smaller and smaller, since heat destabilizes the 
operation of the devices. Therefore proposals for the spintronics devices without creation of charge current flow would 
, be greatly desirable. 

t-H ' Some spin-based devices have been already proposed^. However, one of the most effective way to create a spin 
^sO I current is using semiconductor quantum dots^^iSiZiS. Martin et al. proposed a scheme for electrical detection 
of electron spin resonance(ESR) of an electron trap£. This idea intrigues subsequent studies on how to achieve a 
desirable spin current through a quantum dot or multiple quantum dots either in the strong Coulomb regime^i^ or 
in the Kondo regime^. The basic idea is first Zeeman splitting the dot level by a perpendicular constant magnetic 
field and then pumping electron from the low-lying spin-up state to the high-lying spin-down state by an oscillating 
magnetic field. We notice that special parameters are chosen to facilitate the studies£&££. It is the purpose of this 
work to relax the parameter constraints to gain a general physical picture about spin pump effect in a quantum dot 
Q , by a time-varying magnetic field. 

O ■ The quantum dot spin pump consists of a quantum dot coupled to two electronic leads by tunnel barriers. No 
voltage bias is applied to the device to allow for a direct current through the quantum dot. To realize a spin 
pump with such a system, one apply a perpendicular constant magnetic field Bqz and a lateral time- varying field 
B r f[cos(u> r ft)x + sin(u) r ft)y] to the quantum dot. Due to Zeeman effect in the presence of a constant magnetic field 
Bqz, the dot level will be split into two spin-dependent levels: Sd — ► s a = Sd — (tEz{= qiibBq) /2, where g, fiB are the 
effective g factor and Bohr magneton of the quantum dot. The oscillating magnetic field serves a machine to pump 
low-lying spin-up electrons to the high-lying spin-down level, which is described by a Hamiltonian H r f given below. 
We consider the case that the Coulomb interaction between electrons inside the dot is strong enough to prohibit 
double occupation of the dot. A spin- up electron tunnels into the quantum dot from either the left or right lead, and 
tunnels out of it after being pumped to the high-lying state. This process persists repeatedly and a steady charge 
and spin current is generated. 
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The Hamiltonian of this device is 

a— l 7 r a—l,r 

H = ^^e a d) a d a + ^ £ak<jc' akcr c a k<, + ^ t akcr (& ak(T d a + d\c aka ) + H rf (t), (1) 

u ka ka 

where d^.(d a ) creates (annihilates) an electron with spin a =|,J. in the quantum dot at the level e a) c^ ak<J (c aklJ ) 
creates(annihilates) a spin-cr electron in the lead a = I, r. The third term in Eq. (1) describes tunneling between 
the dot and the leads, while the last term H r f denotes the pumping mechanism and can be written as H r f(t) = 
fl R (<tfdie iu '-t t + d J f d T e-^*)/2, where il R = gfisBrf /2{h is set to be unity throughout this work) is the Rabi 
frequency. 

A gate voltage must be applied to the quantum dot and so adjusted that the chemical potential of the leads lies 
between the spin-up and spin-down levels of the dot. Initially the system lies in its ground state \G) with electrons 
filling up to the chemical potential in the leads and without electron inside the dot. When tunneling is turned on, the 
wave function of the whole system can be written in the following form 

a—l,r a : fi—l : r 

|*(t)> = {b (t)+ Kki(t)d\ 

CafcT + b ak i(t)d^c ak ^] + ^ b a p(t)c^ k , ^c ak ^ + 

a,/3,7— Z,r 

[^/3 7 T(*)4 C 7fc"i Cafc T C /3 fe 'T + b aPll {t)d\c\ klli C ak] Ci3k'\] H }\G) , 

where b(t)'s denote the probability amplitudes for finding the system in the corresponding states at time t, and 
initially all the b(0)'s except bo(0) are zero. 

Now we introduce the reduced density matrix pij spanned in the Fock space of the quantum dot: |0) — 1 empty 
state, | j) — the spin-up state is occupied ,11)—^ the spin-down state is occupied. The diagonal elements of 
the density matrix pn give the probabilities of finding the dot being either empty or occupied by a spin-cr electron, 
while the off-diagonal elements describe coherent superposition of the spin-up and spin-down states. The density 
matrix p^ can be obtained by tracing out the degrees of freedom of the leads in the full density matrix p^ = 
Sn !T m -n r n Tl p -™ !T '™ u '™ T ' T '™ rJ - \ where p^ l T' ni i' nr f' nr ^ represent the probabilities of finding the dot in the state ij 
with ni\ and n r j spin-up electrons tunneling out of the left and right leads, and riiy and n r ^ spin-down electrons 

tunneling into the left and right leads. We find p o = \bo(t)\ 2 + J2k'k~^ T \ b ap(t)\ 2 H , /W = Y,k~ l ' T \b a ka{t)\ 2 + 

Efe,f'',l»'' r \b a 0-r<y{t)\ 2 + ■ ■ ; p n = TT lT b a k]{t)b* aki {t) + J2tfk'T' l ' r b ^l(t)b* al3ll (t) + ■■■■ Solving the schrodinger 
equation i^p 1 = H\V(t)) results in an infinite set of coupled linear differential equations for b(t)'s, which can be 
finally transformed into an infinite set of algebraic equation o 9 ' 10 after the Laplace transform b(E) = J °° dtb(t)e lEt . 
Performing inverse Laplace transform and summing up the relevant terms, we obtain the follow quantum rate equations 
for the density matrix p^ 

Poo 
PTT 

Pll 

PIT 

Pn 

where F QCT = 2irJ2k \t a ka\ 2 S{uJ 

the dot and the leads. We have introduced phenomenologically an additional relaxation term 1/Tj_ to describe the 
transverse spin relaxation process. 

Current l aa is calculated as the evolution rate of electron number tunneling into or out of the lead a: I aa = 
eN ary (t) — n aa p^ n '' n ' l ' nrurirl \ In the stationary state limit, p\j = 0, one finds the following expressions for the 

charge current 7° = 7 Q j + 7 Q j and the spin current 7* = I a i — I a -\ flowing into the left or right lead 



= -(r JT + r rT )poo + (r u + r ri )p u , (2) 

= (r ;T +r rT )poo + *^(e^V Ti -e-^ViT), (3) 

= -(F u + T rl )p u i^( e *-r/* pT1 _ e"^VlT), (4) 

= (iE z ^±^i _L )piT + i^e^\p n - p TT ), (5) 

= {-iE z ^i±En _L )PU + i^ e -i^ t{Pn _ m); (6) 



— £aka) are the line- width functions characterizing the coupling strength between 



if = i c r = en%(r n r rl - r u r rT )e(B z> n fl> w P /,T L ,r aff ) J 



(7) 
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if = -en^(2r, T r u + r !T r ri + r u r rT )e(£7 z ,n ii ,w r/ ,r L> r a<7 ), 
i^ = ef2 i j(2r J .-[-r r | + r r ^rj| + r r ^T^)@(E z , o. R ,u) T f, T_\_,T aa ), 



(8) 
(9) 



where the resonance function is 



Q(Ez,£Ir, u) r f,T±, T aa ) 

T(r ± ,r Qi ) 



T(Tx,r Q .4)/A(E z ,fifl,w r /,Tj_,r Q!0 .), 

a{e z - u;rf) 2 (T n + r rT )(r u + r ri ) + n|(2r n + 
2r rT + r u + r ri )T + (r ZT + r rT )(r u + r ri )r 2 . 



(10) 



(ii) 



A(E Z , Or, u r f,T±, T aa ) 



(12) 



We notice that Dong et ali£ considered a particular case with all line-width functions being the same, and drew the 
conclusions of zero charge current and same values for the spin current in the two leads. We see from Eq. (7) that 
the pumped charge current in the left and right leads has the same magnitude and direction no matter what the 
parameter values. It is a result of current conservation. Charge current disappears when the line-width functions 
satisfy the relation T^/T^ = r r |/r r |. While the condition to have the same magnitude for the spin current in the 
left and right leads is IVr^ = IY-flY^. In the spin-independent tunneling case T a ^ = T a i, charge current is always 
zero, and a maximal spin current can be expected when the time-varying field is resonantly coupled to the dot, i.e., 
Lu r f — Ez, and the dot is coupled to the leads in an extremely asymmetric way. The ratio \If /J*| between the 
magnitudes of spin current in the left and right leads is directly proportional to the coupling asymmetry factor T; /T r 
in the spin-dependent tunneling case. 

In summary, we have derived explicit expressions for the charge and spin current in terms of setup parameters of a 
quantum dot pump device, and discussed a possible way of optimizing the relevant parameters to achieve a maximal 
spin current without charge current flow in the leads. The influences of finite coulomb interaction and spin- flip process 
on the charge and spin current in such a device in the nonequilibrium situation are expected to be more interesting 
and the study is underway. 
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